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Abstract. We present first worldline analytical and numerical results for the 
nontrivial interplay between geometry and temperature dependencies of the Casimir 
effect. We show that the temperature dependence of the Casimir force can be 
significantly larger for open geometries (e.g., perpendicular plates) than for closed 
j5 ' geometries (e.g., parallel plates). For surface separations in the experimentally relevant 

^ , range, the thermal correction for the perpendicular-plates configuration exhibits a 

stronger parameter dependence and exceeds that for parallel plates by an order of 
magnitude at room temperature. This effect can be attributed to the fact that the 
fluctuation spectrum for closed geometries is gapped, inhibiting the thermal excitation 
^ ' of modes at low temperatures. By contrast, open geometries support a thermal 

CO , excitation of the low- lying modes in the gapless spectrum already at low temperatures. 

O 

^ ■ 1. Introduction 

l> 

^ ' The Casimir effect [1] is a paradigm for fiuctuation-induced phenomena. Casimir forces 
between mesoscopic or even macroscopic objects which result from fluctuations of the 
ubiquitous radiation field or of the charge distribution on the objects inspire many 
branches of physics, ranging from mathematical to applied physics, see [2] for reviews. 
Since fluctuations usually occur on all momentum or length scales, they encode both 
local as well as global properties of a given system. In the case of the Casimir 
effect, the resulting force is influenced by localized properties of the involved objects 
such as surface roughness as well as by the global geometry of a given configuration. 
From a technical perspective, localized properties can often be taken into account by 
perturbative methods owing to a separation of scales: e.g., the corrugation wavelength 
and amplitude are usually much smaller than the object's separation distance. But 
global properties such as geometry or curvature dependencies generally require a full 
understanding of the fluctuation spectrum in a given configuration. 

Recent years have witnessed the development of a variety of new field-theoretical 
methods for understanding and computing fluctuation phenomena. So far, only 
phenomenological recipes had early been developed for more complex Casimir 
geometries, such as the proximity force approximation (PFA) [3]. For the special case 
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of Casimir forces between compact objects, field-theoretic results in asymptotic limits 
had been worked out [H [5] . A first field-theoretic study of the experimentally important 
configuration of a sphere above a plate [6] was performed in [7] based on a semiclassical 
expansion. A constrained functional-integral approach, as first introduced in [8] for the 
parallel-plate case, was further developed for corrugated surfaces in [9]. 

The sphere-plate as well as cylinder-plate configuration [10] was also used as a 
first example for the worldline approach to the Casimir effect [11], which is based 
on a mapping of field-theoretic fluctuation averages onto quantum-mechanical path 
integrals. This technique is rooted in the string-inspired approach to quantum field 
theory which is particularly powerful for the computation of amplitudes and effective 
actions in background fields [12]. For arbitrary backgrounds, the path integral over 
the worldlines representing the spacetime trajectories of the quantum fluctuations 
can straightforwardly be computed by Monte Carlo methods, as first demonstrated 
in [13]. The particular advantage of the approach arises from the fact that the 
computational algorithm can be formulated independently of the background. This 
makes the approach so valuable for Casimir problems, where a given surface geometry 
constitutes the background for the fluctuations. The resulting technical simplifications 
become particularly transparent for fluctuations obeying Dirichlet boundary conditions 
(b.c), where high-precision computations have been performed, e.g., for the sphere-plate 
and cylinder-plate case [HI [151 HB] • 

A number of further first-principles approaches for arbitrary Casimir geometries 
have been developed and successfully apphed in recent years. The constraint functional- 
integral approach has been extended to general dispersive forces between deformed 
media [17]. In particular, approaches based on scattering theory have proved most 
successful, starting with an exact study of the sphere-plate configuration with Dirichlet 
b.c. [18]. Scattering theory also lead to a solution for the cylinder-plate case which, 
as a waveguide configuration, allowed for a study of the case with real electromagnetic 
b.c. [19]. These scattering tools have been further developed to facilitate an analytical 
computation of the important small-curvature expansion [20]. For configurations with 
compact objects, new scattering formulations have recently been found which separate 
the problem into the scattering off the single objects on the one hand and a propagation 
of the fluctuation information between the objects on the other hand [2T1 [22] : in 
particular, electromagnetic b.c. for real materials can conveniently be addressed with a 
new formulation which emphasizes the charge fluctuations on the surfaces [22]. Let us 
also mention the combination of scattering theory with a perturbative expansion that 
has recently allowed to study geometry effects beyond the PFA [23]. Scattering theory 
is also a valuable tool for analyzing Casimir self-energies [23]. Finally, direct mode 
summation has also successfully been applied to nontrivial geometries [25] . 

In a real Casimir experiment, further properties such as finite conductivity, surface 
roughness and finite temperature have to be accounted for in addition to the geometry. 
Generically, these corrections do not factorize but reveal a nontrivial interplay. For 
instance, the interplay between dielectric material properties and finite temperature 
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PB] still seems insufficiently understood and has lead to a long-standing controversy 
[271 [28] • In the present work, we confine ourselves to the ideal Casimir effect where 
this controversy does not exists; but even in the ideal limit, the interplay between 
geometry and temperature can be substantial, as demonstrated below. The difference 
is not only of quantitative nature, but arises from the underlying spectral properties 
of the fluctuations, as first pointed out by Jaffe and Scardicchio [29]. In the familiar 
parallel-plate case, the nontrivial part of the spectrum transverse to the plates exhibits 
a gap of wave number /cgap = ir/a, where a is the plate separation. At temperatures 
T smaller than this gap, the relevant fiuctuation modes are hardly excited, implying a 
suppression of the thermal corrections; the leading small-temperature contribution to 
the parallel-plates Casimir force scales like (aT)^. Geometries with a gap in the relevant 
part of the excitation spectrum are called closed. Following the same line of argument, 
we expect a suppression for thermal effects for all closed geometries. 

By contrast, there is no reason for this strong suppression of thermal corrections in 
open geometries which do not have a gap in the fluctuation spectrum. The sphere- 
plate or cylinder-plate cases belong to this class. For open geometries, there are 
always Casimir-relevant modes in the fluctuation spectrum that can be excited at any 
small value of the temperature. Hence, we expect a much stronger dependence on the 
temperature, e.g., (aT)" with < a < 4, and thus a potentially much stronger thermal 
contribution in the experimentally relevant parameter range aT ~ 0.01 . . . 0.1. 

So far, no flrst-principle calculation has been able to conflrm this expectation, since 
generic asymptotic-limit considerations and standard approximations typically break 
down in the relevant parameter range, as already emphasized in [29]. For instance, 
the exact solution for the cylinder-plate case allowed for an explicit temperature study 
of the limit of small cylinder radius, R -C a,/5, where (3 = 1/T. In this limit, a 
log-modifled (aT)^ correction is obtained for the dominant part of the spectrum with 
Dirichlet b.c. [19]. This result suggests that the low- lying thermal excitations with long 
wavelength are not suppressed by a gap but by the smallness of the cylinder radius 
required by the asymptotic-limit considerations. 

Also, the use of recipes such as the PFA can lead to a different scaling, such as 
a (aT)^ law for the sphere-plate case [6l |30]. Whereas the PFA at zero temperature 
is justiflable in the low-curvature limit, a <C i?, [3, [HI [31], [15], PFA-deduced thermal 
corrections can be problematic: at small temperatures with thermal wavelength much 
larger than the minimal surface separation aT -C 1, the thermal excitations can be more 
sensitive to the curvature radius than the vacuum fluctuations. Even worse, the PFA 
uses the parallel-plate formula, and hence a gapped spectrum, as an input and thus 
misses the important difference arising from an open geometry. 

In this work, we present flrst evidence for a strong thermal correction to a Casimir 
force law for an open geometry using worldline numerics. As a paradigmatic example, 
we use the conflguration of a semi-inflnite half-plate perpendicularly above an inflnite 
plate (cf. Fig. [1]), imposing Dirichlet b.c. for the fluctuations of a real scalar fleld. This 
configuration belongs to a set of cases, revealing a universal force law determined by 
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Figure 1. Left panel: sketch of the paraUel-plate configuration (taken from [32]). 
Right panel: sketch of the finite-temperature spacetime; a worldlinc can wind around 
the compactified time dimension. 



dimensionality, which has first been investigated in the context of Casimir edge effects 
[32]. Since the configuration has only one length scale which is the distance a between 
the edge of the half-plate and the infinite plate, the interplay of the gapless fluctuation 
spectrum with thermal excitations is not disturbed by other length scales, resulting in a 
clean thermal signature of an open geometry. Our worldline calculations yield a thermal 
correction obeying an [aT)^ force law at low temperature. This implies a substantial 
increase of the thermal contribution compared to those for a closed geometry. 

The fact that geometry and temperature exhibit such a nontrivial interplay in 
Casimir systems, resulting in "geothermal" Casimir phenomena[j|, is another peculiar 
feature that should be added to the long list of peculiarities of the Casimir effect; it 
clearly deserves further investigation. 

2. Worldline approach to the Casimir effect at finite temperature 

Let us briefly summarize the worldline approach to the Casimir effect. More detailed 
descriptions and derivations from first principles can be found in flTi [T6] . We consider 
the Casimir interaction energy, serving as a potential energy for the force, for two rigid 
objects with surfaces Si and S2. For a massless scalar field with Dirichlet boundaries 
in D = 3 + 1, the worldline representation of the Casimir interaction energy is given by 



Here, the worldline functional Os[x(r)] = 1 if the path x(r) intersects the surface 
S = Si U S2 in both parts Si and S2, and OE[x(r)] = otherwise. 

This compact formula has an intuitive interpretation: the worldlines can be viewed 
as the spacetime trajectories of the quantum fluctuations of the scalar field. Any 
worldline that intersects the surfaces does not satisfy Dirichlet boundary conditions. 

I We introduce the attribute "geothermal" here, since it directly describes the source and nature of 
this phenomenon. No link exists between the physics discussed here and, e.g., geothermal heat-pumps 
etc. dealt with in the geological sciences. 




(1) 
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All worldlines that intersect both surfaces thus should be removed from the ensemble 
of allowed fluctuations, thereby contributing to the negative Casimir interaction energy. 
The auxiliary integration parameter T, the so-called propertime, effectively governs the 
extent of a worldline in spacetime. Large T correspond to IR fluctuations with large 
worldlines, small T to UV fluctuations. 

The expectation value in Eq. ([1]) has to be taken with respect to the ensemble of 
worldlines that obeys a Gaufiian velocity distribution 

{...)=[ ... e-^^o^ /" Dxe-3/o^'^-*'W, (2) 

where the worldlines have a common center of mass xqm- At zero temperature, the time 
component of the worldlines cancels out for static objects, hence the straightforward 
Monte Carlo computation of Eqs. ([I]) and ([2]) can be restricted to the spatial part. 

Finite temperature can now easily be implemented with the aid of the Matsubara 
formalism, and also the technical changes of the numerical algorithm are only minor: 
The Euclidean time, say along the -Dth direction, is compactified to the interval [0, /3] 
with periodic boundary conditions for bosonic fluctuations. As a consequence, the 
worldlines can also wind around the time dimension, see Fig. [H It is convenient to write 
a given loop x(r) with winding number n as sum of a loop with no winding, x(r), and 
a translation in time running from zero to n(3 with constant speed, 

T 

x^{t) = x^(r) + np—S^D- (3) 

The path integral over the different winding number sectors labeled by n factorizes for 
static configurations, yielding 

Pa;e-/o"'^-4 ■■■ = E ^""^ f Vxe-fo^^^ .... (4) 

x{o)=x{T) Jx(o)=H'r) 

The worldline representation of the Casimir interaction energy for the Dirichlet scalar 
at finite temperature thus reads 

1 1 dT / °^ 2 2\ r 

^Casimir = "2?^ / T^\Y. ^"^^j / ^'^^M (es[x(r)]) . (5) 
^ ' n=-oo 

Whereas the worldline expectation value remains identical to the one at zero 
temperature, the winding sum re-weights the propertime integrand: larger temperature 
emphasizes smaller propertimes and vice versa. This confirms the expectation that 
thermal corrections at low temperature are dominated by long wavelength fluctuations 
which in our case correspond to worldlines with a large spatial extent. 

It is important to note that -Ecasimir is normalized such that -Ecasimir — > for infinite 
distances a oo. Hence, -Ecasimir can differ from the thermodynamic free energy by 
thermal corrections to the self-energies of the single surfaces. The latter is distance- 
independent and thus does not contribute to the Casimir force. 
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Figure 2. Temperature dependence of the Casimir force between two parallel plates. 
The ratio of the force at temperature T and at zero temperature is plotted versus the 
temperature in units of the plate distance a. The rectangle at low temperature marks 
the region magnified in Fig. [3l For the worldline numerical result, we have employed 
800 loops each with 1 000 000 points. 



2.1. Parallel plates 

As a test, let us consider the two parallel plates separated by a distance a along the 
z axis. Interchanging expectation value and zqm integration in Eq. Q, we encounter 

oo 

dzcMQE[^] = {Vfl-a)e{^/Tl-a), (6) 

■oo 

where / denotes the dimensionless extent of the given worldline in z direction measured 
in units of a/T ; cf. [16]. Differentiating Eq. ([5j) by —8/ da yields the Casimir force 

where A is the (infinite) area of the plates. Figure [2] shows the numerical result for the 
Casimir interaction energy ([7]), corresponding to the distant-dependent part of the free 
energy, normalized to the zero-temperature result. For comparison, the analytic result. 



7 - 
Casimir r> o 

2 oa 



— 2~ Z-^ 7 ^3 (coth(27rmaT) + 27rmaT csch^(27rmaT)) 

m=l ^ ' 



(8) 



is also shown, see e.g. |33)§l . Both results agree satisfactorily. 

Incidentally, the leading thermal correction can be obtained analytically from the 
worldline representation (JTj): for (aT)^ <^ 1 (and n 7^ 0), the propertime integrand is 
dominated by large T, hence the lower bound can safely be set to zero. This results in 
the well-known leading thermal correction AF(T) = — (7r^/90)AT^, which can also be 
understood as an excluded volume effect: thermal modes are excluded from the region 
between the two plates which thus does not contribute to the Stefan-Boltzmann law. 

§ We use half the value of [33] which is derived for the electromagnetic field with two degrees of freedom. 
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Figure 3. Temperature dependence of the Casimir force for two perpendicular plates 
compared to the parallel-plates result. The ratio of the force at temperature T and at 
zero temperature is plotted versus the temperature in units of the plate distance a. The 
plot shows the small-temperature range of Fig. |2l At experimentally relevant large- 
separation values of To (vertical line), the temperature correction for the perpendicular 
plates is ~ 6%, which should be compared with ~ 0.7% for the parallel plates. For 
the worldline numeric results we have employed 800 loops each with 1 000 000 points. 
The error bars represent the statistical error. 



2.2. Perpendicular plates 

We now consider the perpendicular-plate configuration introduced above. Again, we 
can perform the Zcu integration first, yielding for the force 

^0-.,, - ^ (/ X,,^,, ^ ( f ^ -) ) . (3) 

where L is the (infinite) length of the system along the edge. Here, /(^) denotes the 
dimensionless extent of the given worldline in z direction as seen by the configuration in 
units of a/T ; it depends on the position ^ of the worldline normal to the perpendicular 
plate which is also measured in units of VT ; for details, see [M]. Figure [3] compares 
the resulting temperature correction with that of the parallel-plates case in the small 
temperature range of Fig. [2l In contrast to the weak (aT)^ dependence of the parallel- 
plates result, the Casimir interaction energy for the perpendicular plates shows a strong 
increase with temperature. For typical experimental values at larger distance a = 1.5/im 
and room temperature, the temperature correction is about 6%. At the same distance 
and temperature, the temperature effect for the parallel plates is 0.7%. The open 
geometry therefore exhibits a thermal correction which is an order of magnitude larger 
than the closed parallel-plates case. 

The leading thermal correction can again be computed analytically from the 
worldline representation by extending the lower bound of the T integral to and 
using (J d^) = (l) = y/n. Here, / denotes the extension of the loop perpendicular to the 
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semi-infinite plate We obtain 

C(3) L 

Fcasimir(T) ^ Fcasimir,T=0 " ^ ^ i^Tf, for (flT) < 1, (10) 

which is confirmed by the full numerical result over the whole range of temperatures 
shown in Fig. [3l We note that the low-temperature scaling of the thermal correction for 
this open geometry cannot be understood as an excluded- volume effect. 

3. Conclusions 

We have presented analytical as well as numerical results for the nontrivial interplay 
between geometry and finite temperature for the Casimir effect in open geometries. For 
the first time, we have shown that the gapless nature of the fluctuation spectrum leads 
to a strong enhancement of the thermal correction to the Casimir force. Our numerical 
data for the perpendicular-plate case with Dirichlet b.c. confirms our analytically derived 
(aT)^ force law at low temperatures. This should be compared to the weaker (aT)^ 
dependence of the parallel-plates case where a gap in the spectrum suppresses the 
thermal correction. This calls urgently for further first-principles computations for other 
open geometries such as the sphere-plate case. 

It is a pleasure to thank Michael Bordag and his team for the organization of the 
QFEXT07 workshop and for creating such a stimulating atmosphere. This work was 
supported by the DFG under Gi 328/1-3 (Emmy-Noether program) and Gi 328/3-2. 
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